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I. INTRODUCTION 


During last years string four-dimensional dilatonic black holes attracted much attention. 
As this type of black holes is the solution of the string theory in its low energy limit, therefore, 
by studying these solutions one can hope to clarify some important unsolved problems of 
modern theoretical physics, for example, what is the endpoint of the black hole evaporation, 
the quantum coherence and black hole thermodynamics problems M. SO on. After the 
appearance of the Gibbons-Maeda-Garhncle-Horowits-Strominger (GM-GHS) solution a 
great interest to the investigation of the higher order curvature corrections in the Einstein- 
dilaton (Yang-Mills) lagrangian arised |]5|-|Tll. It is now very actual problem because in the 
regions where the curvature of the space-time increases the role of the a' expansion terms 


grows. The general form of this expansion is not investigated well yet |]^ . So as a hrst step 
researches study the contribution of the second order curvature corrections. It was found 
that the black hole solutions with the {a'Y terms generalize the wellknown Schwarzschild 
solution and the GM-GHS one. The problem is to find new feathers of the black hole 
solutions which were introduced by the string theory. 

Some researchers study only the simplified bosonic part of the low energy string action 
taken in the form 
1 


^ = 


(TXs/^ 


Ibvr 

- + Xe-'^^ScB 


mpi(—R + 


( 1 ) 


where i? is a Ricci scalar; 0 is a dilaton held; mpi is the Planck mass; F^^, is the Maxwell 
held; A = a 'is the string coupling parameter describing the Gauss-Bonnet (GB) term 
contribution {Sgb = —ARijR^^ + R"^) to the action (|l|). It was found that the black 

hole solution of such action (with or without the Maxwell term) does exist and provides 
the non-trivial dilatonic hair [p|-[T0[|. Further, it was established that the modihcation of the 
solution by the second order curvature corrections became non-vanishing only if the black 
hole size (mass) was small enough (the case of the large value of the coupling constant A). 
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It is necessary to emphasis that the differential equations in such class models have a very 
complicated structure. So, there is no possibility for the direct analytical solving them. 
Hence, the researchers have to use perturbative |]T0| , p!^ , p!5| or numerical methods. 

In our previous papers [0,|| black hole type solutions of the action (|^) were found by the 
special numerical method from infinity up to a some particular point inside the event horizon 
Th- It was shown that in the case when the magnetic charge q was rather small (or vanished) 
a singular “tube” (in t direction) with the topology S‘^ x has appeared inside the black 
hole. Asymptotically flat solution occurs from infinity up to this singular “tube” with the 
radius r^. It takes place in the range of the magnetic charge g to be 0 < g < gcr, where qcr 
is a new critical charge value appearing in second order curvature gravity |§]. In the case 
of g to be large enough (g^^ < g < M\/2 in GHS gauge, see @]) the solution occurs up to 
zero point where the dilatonic function 0 diverges as in the GM-GHS case. In this paper 
we are going to show that all the particular points (and their main feathers) of the black 
hole solution obtained from the set of implicit ordinary differential equations are defined by 
the various kinds of zeros of the main system determinant. In addition some new feathers 
of these solutions are discussed. 

The structure of the paper is the following. Sec. II deals with the main equations, in Sec. 
Ill we briefly remind our previous results and discuss them in the light of the conclusions of 
T^, in Sec. IV some new physical feathers of the singularity are presented and Sec. V 


contains the main conclusions. 


II. FIELD EQUATIONS 

Our purpose is to study physical and mathematical feathers of the black hole solutions 
(static, spherically symmetric case). Therefore, the most convenient choice of metric is 

= A{r)dt^ — ~ + sin^ Odcp^), 

where functions A, a and / depend only on a radial coordinate r. Various types of the 
metric gauges were used in our investigations. The most convenient of them are, so-called. 
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“curvature gauge” (/ = r, used in 0) and GHS gauge (a = 1, used in 0). Corresponding 
field (Einstein) equations can be rewritten in the matrix form 


an/S" + a^f + 0*30" = bi, GHS gauge 

OjiA" + di 2 (j' + 0*30" = bi, “curvature gauge” (2) 

where i = 1,2,3 and the coefficients o*j and bi are equal to (in the GHS gauge) 

Oil = 0, 

oi2 = + 4e“^'^A0' 2A/', 

013 = de-^-^A (A/'^-l), 

021 = 'OTp;/ + 4e“^'^A0' 2A/', 

022 = mil 2A + 4e-2'^A0' 2AA', 

023 = de-^-^A 2AA7', 

031 = 4e-2<^A (Af'-l), 

032 = de-^-^A 2AA7', 

033 = (-2)m|;A/^ 

h = mhn^f 

+ de-^-^A (Af ^ - 1) 2(0')^ 

62 = (-2)m|,(A7' + A/(0')7 
+ de'^'^A 2AA7' 2(0')^ 

+ 2e-2V(l/f) 

- de-^-^Af 2(A')V', 

63 = 2m?,77 (A7 + 2A/') 

+ 2e-2 V(l/f) 

- de-^'^A (A')2/'^- 
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The remaining SS/6f = 0 additional equation fulfills automatically elsewhere on the solution 
trajectory. 

The Arnowitt-Diezer-Misner (ADM) mass M and the dilaton charge D are dehned by 
the asymptotic expansions for A, / and (j) 


A 

f 

0 


2M 

r 




r 



(3) 


III. INVESTIGATION OF THE PARTICULAR POINTS 


For integrating outside and inside the event horizon, a method based on integrating over 
an additional parameter was used and described in detail in [^]. Here we briefly remind the 
main results and discuss their new feathers. 

The results of our numerical integration are depicted in Fig.l. Shown in the Figure 
presents the dependence of metric function A against radial coordinate r with the different 
meanings of the event horizon value Vh (in the curvature gauge). It was found in [Q that the 
solution with the singular turning point Vg exists only if 0 < g < gcr where qcr is the critical 
magnitude of the magnetic charge. The appearance of Qcr is created with the second order 
curvature corrections. While q increases and reaches the qcr value, the solution behavior 
changes and it looks like completely to the GM-GHS one. As one can see from Fig.l, the 
behavior of the above-mentioned metric function A outside the regular horizon has the usual 
form and is similar to the standard Schwarzschild solution (see, for example, |p|,[l0|). Inside 
the horizon the solution exists only down to the value r = r^. Another solution branch 
(additional branch) begins from the value and exists only up to the “singular” horizon r^- 

To check the particular points of the system (0) it is necessary to consider its main 
determinant (here, for simplicity, we use “curvature gauge” and vanishing value of magnetic 
charge; in the GHS gauge with 0 < g < gcr the results are completely the same but the 
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formulas have larger length) 


D 


main 


A 


+ BA+ C 


(4) 


where 


A = (-32)e-^‘^a^A^ 


A 2 i./2 2 2 

4(J 0 mpiT — 


- mil+ 12e-^^A'^'X 

B = (—32)e“^'^cT^A a‘^(f)'mlir^ + 

+ 2e-^^a^Xmli - Se'^'^AVA^ 
C = 32e-^‘^a^X^mli - 2a®m5;r^ + 


+ 64e-^'^a® A'A^mJ^r + 128e-®'^a® A'0'A^ 


Eqs. represent the system of ordinary differential equations in a non-evident form. 
That is why the system can have the peculiarities when its main determinant Dmain 
vanishes. The structure of those peculiarities depends upon the behavior of the eqs. 
near the singular surface Dmain = 0 in phase-space [|^. In the model with asymptotically 
flat solutions only three types of the Dmain degeneracy arise. They are 


(a) : A = 0, C 7^ 0 (“intersection”), 

(b) : AA^ + BA + C = A 7^0, C^O (“turning point”), 

(c) : A = 0, C = 0 (“complicated singularity”), (5) 

The case (j^a) is realized at the regular horizon (see Fig.l) and occurs also in the 
Schwarzschild solution. The asymptotic form of the solutions near the point Vh has the form 

0 

A = di(r - Th) + d 2 {r - Vhf + o((r - 
a = So +Si{r - Th) + o{{r - Th)), 

e"^"^ = 0o(l - 2 * 0i(r - Th) + 2(01 - 02)(r - rhf) + o((r - Thf), (6) 
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where {r — rh) -C 1. Substituting the formulas (|[) into the equations (||), one obtains the fol¬ 
lowing right relations between the expansion coefficients (sq, 0 o and are free independent 
parameters, here we do correct an unfortunate misprint in our work |]^) 

di(^zidi Z2di ^ 3 ) = 0 , ( 7 ) 


where: 

= 24AVo, 

4 2 4 

^3 = rUpir^SQ, 


and the parameter 0 i for di 7 ^ 0 is equal to: 0 i = [(mp;)/(4A(ii0o)] * [r^di — Sq]. 

When di = 0, the metric function A has the double (or higher order) zero. In such a situ¬ 
ation the equation for ^2 {d^, d^^, ...) is a linear algebraic one and there are no asymptotically 
flat branches. 

When di 7 ^ 0 the solution of the black hole type takes place only if the discriminant 
of the equation (]^) is greater or equal to zero and, hence, r\ > 2A0o\/l4 -|- 8 a / 3 . One or 
two branches occurs and always one of them is asymptotically flat. With the supposition 
of 0 OO = 0 (and, as we tested, in this case 1.0 < 0 o < 2 . 0 ) the infinum value of the event 
horizon is 




= V\ 


\ 


4^ 


m 


2 ’ 

PI 


( 8 ) 


The analogous formula in the other interpretation was studied by Kanti et. al. in |^. 

The case (^d) is realized inside the regular horizon at r = Vg < Vh (see Fig.l) as a 
consequence of the intensihcation of the GB-term influence. The solution behavior strongly 
differs from the Schwarzschild one or the GM-GHS one. The D^ain degeneracy of (jpj) 
type reduces to the violation of the uniqueness of the solution at the point r^. Similar 
situations are typical for the systems of the type (|^) in the neighborhood of the singular 
surface D^ain = 0 [M. The asymptotic behavior of both solution branches near the position 
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Vs can be described by the following formulas using the smooth function a as an independent 
variable 

A = ds + d 2 {a - + o{{a - asf), 

r = Vs + r 2 {(T - asf + o{{a - asf), 

exp(-20) = 0^(1 - 2 / 2 ( 0 -- + o((cr - (9) 

where o' — o-^ 1. Free independent parameters are: as, 4>s, fs- After the substitution of 

these expansions to the system (H), one can obtain that ^2 = /2 and the equation for (^ 2 /^ 2 ) 
has the form 

^42/^ + + Zti = 0, (10) 

where y = {d 2 /r 2 ) and the other coefficients are 

Z 4 = mpia‘^dsr‘1 + 4A0s(o-^ - 3ds), 

Z5 = -mpia^Ts, 

Z6 = m%a^(a^-ds). 

Eq.(pi]D may have either no solutions or one or two solutions dependently upon its dis¬ 
criminant magnitude. The case of the unique solution corresponds to the minimal value 
of ’"hmin (eq- (§), see Fig.l curve (c)). The situation with the positive discriminant cor¬ 
responds to the turning point of the solution. Here it is necessary to note that if one 
rewrites the expansions (P) against the expansion parameter (r — r*) -C 1, the result will 
have the following form A = -|- y(r — Tg) + ■ ■ ■ and so on. Hence, there are only two 

branches (because of two possible values of y) can exist near the position r^. They are: 
the asymptotically flat one and the {rsTx) one. No any other solution branches are in the 
neighborhood of rg. The “curvature invariant” RijkiR'^^^ near is equal to (r — r* 1): 

= Ti/(r - rg) + o(l/(r - r,)) —>■ 00 , where Ti = const. The components Tq and 
Tf of the stress-energy tensor Tjf also diverges (as l/\/r — Vg) near the position r*. 

The case (^c) is realized on the singular horizon r^ of the additional branch (r^rj.) (see 
Fig.l, curves (a) and (b)). This branch is not asymptotically flat and, hence, is non-physical. 



The asymptotic form of this solution near the position Tx is shown in [0. “Curvature 
invariant” RijkiR'’^^^ also diverges near the position r^- 

Here it is importantly to stress that the distance between the points Tx and Vh decreases 
while decreasing rh (see Fig. 1). In the limit point dehned by the restriction (^) all particular 
points pour together rh=rs=rx=rhmin and the internal structure of the black hole disappears. 
The case (|]c) is realized in this point and, therefore, the “curvature invariant” diverges. 
Hence, the point Thmin represents the event horizon and the singularity in the same point. 
So, such situation contradicts with the “cosmic censorship” hypothesis ||l7|,[l^ but there 
is a question about its stability. The possibility of the realization of such situation is still 
open. Returning to the formula (^, one should remember that A is the combination of 
the fundamental string constants. That is why this formula can be reinterpreted as the 
restriction to the minimal black hole size (mass) in the given model. This restriction appears 
in the second order curvature gravity and is absent in the minimal Einstein-Schwarzschild 
gravity. This fact can throw an additional light on the problems of the black holes in our 
Universe. 


IV. Rs SINGULARITY 


It is possible to obtain the approximate relation between r^, and A. Substituting the 
Schwarzschild values of the metric functions and the vanishing value of the dilaton charge 
D to the formula (®, one obtains 


r, = A 


Ay/3rh(ps 


m 


( 11 ) 


PI 


Fig. 2 shows the dependencies of the value against the coupling parameter A given by the 
formula (pd]) and by the numerical integration. From the eq. ([TI|) it is possible to conclude 
that the pure Schwarzschild solution is the limit case of our one with = 0. In the case with 
rather small value of A this formula gives the good agreement with the results calculated by 
the numerical integration. While increasing A, the absolute error increases as a consequence 


9 







of ignoring the non-vanishing valnes (1 — a), 0' and so on. It is necessary to point ont that 
the eq. o represents the dependence = const which we snppose to be right because 
after the appropriate selection this constant by hands the agreement between numerical data 
and this formula improves. Eq. (^) shows also that when the influence of the GB term (or 
black hole mass) increases, also increases. 

Further, it is possible to hnd the approximate relation between the qcr, A and M. One 
can rewrite the Dmain in the GHS coordinates and substitute there the GM-GHS values of 
the metric functions A, / and the dilaton function 0 in the following form 


A = 

/ = 
e-20 ^ 

Hence, fAmain 


2M 


t '^2 


M’ 


1 - 


1 

M r 


takes the form (we suppose mp; = 


1 for simplicity) 


D 


main 


T 

r^^M\rM - 


( 12 ) 


where T = T{M, A, r) is the polynomial of the sixteenth order against r (we do not write 
it because of its dangerous length). If one supposes the charge q to vanish in the eq. ( p!2[) 
(hence, the denominator of eq. ([T^ never vanishes) he obtains the formula ([Tl| ) in the form 
rl = SVSXM. 


The denominator in the eq. ([T^) can vanish only in the case of r = /M. If this 

situation happens, Dmain diverges. Gonsequently, all the senior derivatives vanishes and 
then in the point which is also particular one has only the local minimum (maximum). 
So, the condition for the turning point existence is the following: must be situated righter 
relatively to r^, i.e. > r^. The limit condition = r*, is just one for qcr because when 
Tfc > Ts one must have a local minimum (maximum) righter r^. That is why equating each 
other the expressions of the and one obtains the approximate formula for qcr 
/ \ 1/6 

ger = . (13) 
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The last formula analogously to the (pT]) gives the good agreement with the data obtained 
from the numerical calculations only in the case of rather small A. While increasing A 
the absolute error increases as well. This fact can be seen from the Figure 3 showing the 
dependence of the value qcr against the coupling constant A. Here it is necessary to point 
out once more that the eq. represents the right dependency qcr = const A^/®. 

Using the expansions (P), it is possible to study the behavior of the radial time-like and 
isotropic geodesics near the singularity (“curvature gauge” case with vanishing q). Based 


on these expansions and on the technics of the geodesical curve calculations from ref. [|l9 


after integrating over the radial coordinate x = r — -C 1 one obtains the expressions 
for the proper time r(x) and the coordinate time t{x) for the radial time-like and radial 
isotropic geodesics 


r(x) 

= iCiX -h C*2 -|- . . . , 

(time-like) 


t(x) 

= lhC*3X -|- C*4 , 

(time-like) 


r(x) 

= iC*5X -f Ce -l-..., 

(isotropic) 


t(x) 

= iCyX -|- C*8 -h . . . , 

(isotropic) 

(14) 


where Ci are the constant values. 

Comparing the values t{x) and r(x) with the Schwarzschild ones near the regular event 
horizon (see, for example, [0), one can conclude that the behavior of the radial geodesics 
near the singularity differs from the radial geodesics behavior of the Schwarzschild solution 
near the regular horizon. In our model all the above-mentioned values limited, in the 
Schwarzschild model t{x) everywhere diverges. The divergence of the also differs 

in the models considered: 1/x® in the Schwarzschild case near the origin and limited at the 
horizon and only 1/x in our case (x = r — r*). Therefore, one can suppose that the r^- 
singularity is a week one (according to the Clarke classification |^) and it can be removed 
by the appropriate extension of the metric. According to the Propositions 8.2.2 and 8.2.3 


11 





from [201, this procedure is not forbidden because near the position v = Vg “criterion’" 


V v' 

K{y) = J dv' J dv"RQo{v") (15) 

0 0 

is limited. So, the question on the possibility of removing the singularity by the appro¬ 
priate metric extension remains open. It is necessary to emphasis that the singular turning 
point Tg appears in various kinds of metric parametrizations as we tested. 


V. CONCLUSIONS 

The black hole solutions generated by the bosonic part of the four-dimensional low energy 
superstring effective action with the second order curvature corrections are discussed in this 
paper. They are obtained by using the special numerical method described in 0. It is 
demonstrated that all the particular points of the solution, namely, regular horizon Vh, 
singular horizon and r^-singularity, are dehned by the various types of zeros of the main 
determinant Dmain, namely, “intersection”, “complicated singularuty” and “turning point” 
corresp ondingly. 

The restriction for a minimal black hole size (mass) is obtained in the frame of the model 
with the vanishing Maxwell held contribution (Einstein-dilaton-Gauss-Bonnet model). This 
minimal black hole size (mass) depends upon the combination of the string fundamental 
constants. 

The approximate formulas for the Vg = rg{X,rh) at vanishing q and Qcr = qcr{X,M) are 
found. These formulas have a perturbative nature and shows that the effects of second order 
curvature corrections become non-vanishing at rather small sizes of black holes. 

The behavior of the radial time-like and radial isotropic geodesics is studied. According 
to this and other criterions Vg singularity is not “strong”. We have no arguments that the 
Ts-singularity is coordinate one (and can be removed by the appropriate metric extension), 
but nothing forbids this. This question is still open. 
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FIGURES 
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FIG. 1. The dependence of the metric function A versus the radial coordinate r at the 
different values of the event horizon ru- The “curvature gauge” with the vanishing contribution 
of the magnetic field is used during the calculations of the data shown in Figure. The curve (a) 
represents the case where is rather large and is equal to 30.0 Plank unit values (P.u.v.). The 
curve (b) shows the changes in the behavior of A(r) when is equal to 7.5 P.u.v. The curve (c) 
represents the boundary case with Vh = J’/imm where all the particular points, namely, rh (regular 
horizon), Vg (singular turning point) and Vx (singular horizon), pour together and the internal 
structure disappears. The curve (d) shows the case where 2M <C rhmin (2M=1.5 P.u.v.) and any 
horizon is absent. The envelope curve (e) shows the position of the points Vg against the different 
values of r^. 
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•A-n*c3 'SH 

FIG. 2. The dependence of the position of the singular turning point r* versus A. Squares 
represent the values of r* calculated from the numerical integration. The solid curve is obtained 
by using formula (11) with = 100.0 P.u.v. and cj)s = 1-2. 
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FIG. 3. The dependence of the position of the critical value of the magnetic charge Qcr versus 
A. Squares represent the values of Qcr calculated from the numerical integration. The solid curve 
is obtained by using formula (13) with M = 5. P.u.v. 
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